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Abstract 

A new mathematical framework is formulated to derive the effective equations of motion for 
the constrained quantum system which possesses an internal clock. In the realm close to classi- 
cal behavior, the quantum evolution is approximated by a finite system of coupled but ordinary 
differential equations adhered to the weakly imposed Hamiltonian constraint. For the simplified 
version of loop quantum cosmology in the Bianchi I model with a free massless scalar filed, the 
resulting effective equations of motion affirm the bouncing scenario predicted by the previous stud- 
ies: The big bang singularity is resolved and replaced by the big bounces, which take place up 
to three times, once in each diagonal direction, whenever the directional density approaches the 
critical value in the regime of Planckian density. It is also revealed that back-reaction arises from 
the quantum corrections and modifies the precise value of the directional density at the bouncing 
epoch. Additionally, as an example of symmetry reduction, we study isotropy emerging from the 
anisotropic Bianchi I model in the context of effective equations of motion. 
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I. INTRODUCTION 

The comprehensive formulation for loop quantum cosmology (LQC) in the k = 
Friedmann-Robertson- Walker (FRW) model (i.e. spatially flat and isotropic) with a free 
massless scalar field has been rigourously investigated [l|, 0], revealing that the big bang 
singularity is resolved and replaced by the big bounce. Later in [3|, based on the same prin- 
ciples but implemented with a more sophisticated quantization scheme, a new Hamiltonian 
constraint was constructed and the big bounce is shown to take place precisely when the 
matter density enters the Planck regime. These attractive results affirm the assertion that 
the quantum geometry of loop quantum gravity (LQG) holds a key to avert the breakdown 
of the classical general relativity and since then have inspired a lot of research on the sim- 
plified (but exactly solvable) model j3, [1] and other extensions such as the k = ±1 FRW 
models [1,0]. To further extend the domain of validity, the framework of [l|, 0, [3| was also 
formulated for the Bianchi I model to include anisotropy [8j. 
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In the Bianchi I model, the similar prediction for the occurrence of big bounces has been 
anticipated in the analytic investigation of jl] and later verified in the numerical computation 
of 0] (for a slightly simplified version of js}). Meanwhile, much effort has been made to study 



the Bianchi I LQC at the level of effective dynamics 

Hal 

111 . Il2| . Not only do the results of 



the effective dynamics agree with those in [s| and 0] but more intuitive pictures are obtained 
in the semiclassical approach. Thus, the study of effective dynamics can give us valuable 
insights into how and why the big bounces take place, even if the definite conclusion for the 
occurrence of big bounces has been drawn directly from the quantum theory. 

However, the validity of the effective dynamics adopted in [Io|, [HI, 12| is only heuristic 
and remains to be justified. A systematic methodology to investigate the effective equations 
of motion for semiclassical states and the back-reaction resulting from quantum effects is on 
demand. 

A mathematical framework has been developed in 13j to study the evolution of quantum 
systems in the realm close to classical behavior. In this framework, the behavior of a wave- 
function subject to a partial differential equation (Schrodinger equation) is approximated by 
finitely many variables subject to a system of coupled but ordinary differential equations, 
giving a technically easier and more intuitive formulation to extract physical information 
of the semiclassical behavior amended by the back-reaction. This formalism has been used 
to address a variety of issues for quantum gravity and quantum cosmology [3, 15]. In 



particular, the bouncing scenario of LQC in the k = FRW model coupled to a free 



massless scalar has been analyzed from this perspective [16j. (An alternativ e ap proach for 
the effective equations of LQC in the k = FRW model was also studied in [171] .) 

In 16], to apply the formulation of [l3[ to LQC in the k = FRW model, the constrained 
quantum system is viewed as an ordinary (unconstrained) quantum system by treating the 
scalar field <ft as the time variable in the very beginning (referred to as "Approach I" in 
this paper). Because of some technical difficulties, Approach I is inept for the Bianchi I 
model. In this paper, we propose a second approach (referred to as "Approach II") as a 
new framework to work with the effective equations of motion for the constrained quantum 
system that possesses an internal clock. In Approach II, the constrained quantum system 
is first treated as if there was no constraint and the formulation of [13j for unconstrained 
systems is adopted to obtain the effective equations of motion with respect to the coordinate 
time without directly treating (f> as the time variable. On top of the finite system of effective 
equations of motion, the Hamiltonian constraint is then weakly imposed up to a certain order 
to further relate the finite set of variables. The technical and conceptual difficulties of the 
standard procedure to quantize a Hamiltonian constrained system are avoided in Approach 
II, the formulation of which can be intuitively understood as the natural extension of the 
classical Hamilton's equations subject to the classical Hamiltonian constraint. Additionally, 
this new approach has an extra merit that measurement of the internal time <fi is on the 
equal footing as that of other observables (namely, the notions of uncertainty and so on for 
are retained), which is an attractive feature for the relativistic quantum mechanics that 
is still missing in the standard quantum theory for constrained systems. With the technical 
and conceptual virtues, Approach II per se could be used as a sound quantum theory for 
the constrained system if an internal time can be identified. 

With the new formulation at hand, we investigate the effective equations of motion for 
the Wheeler-DeWitt (WDW) theory in the Bianchi I model with a free massless scalar. The 
results of Approach II agree with those obtained in the fully developed quantum theory 
in js| for the case without off-diagonal squeezing, thus evidencing the viability of the new 
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formalism. Meanwhile, it is also revealed that back-reaction arises from the off-diagonal 
squeezing. 

The next step is to apply Approach II to a simplified version of LQC in the Bianchi 
model. The resulting effective equations of motion affirm the bouncing scenario anticipated 
in [s| and predicted in [To], 12]: The big bang is resolved and replaced by the big bounces, 
which take place up to three times, once in each diagonal direction, whenever the directional 
density qj approaches the critical value £>/ >C rit m the regime of Planckian density, but the 
precise value of £?/ iCr it is modified by the quantum corrections. Moreover, the off-diagonal 
squeezing gives rise to further back-reaction, which couples the evolution in three diagonal 
directions. 

Additionally, Approach II offers a language to describe the symmetry reduction in terms 
of effective equations of motion. In particular, it is notable that the anisotropic Bianchi I 
model admits the effective solution which exhibits isotropy, but such an effective solution is 
slightly different from that directly obtained in the isotropic model. This observation may 
help us understand more about the issue of symmetry reduction. 

This paper is organized as follows. In Sec. [TTl we first briefly review the mathematical 
framework developed in [l3| and then introduce Approach I and formulate Approach II. In 
Sec. [In] and Sec. HV] respectively, Approach II is adopted to investigate the WDW theory 
and simplified LQC in the Bianchi I model. 1 The issue of symmetry reduction is studied in 
Sec. [Vj Finally, the results are summarized and discussed in Sec. |Vl] 



II. EFFECTIVE THEORY OF QUANTUM SYSTEMS 

In this section, we first briefly review the effective theory developed in [l3| and then 
present two modified strategies — Approach I and Approach II — to study the constrained 
quantum system with an internal clock. Approach II is devised as a new formulation and 
its viability and virtues are discussed. 



A. Effective equations for ordinary quantum systems 



The effective theory of [13] is used to approximate the evolution of a wave-function subject 
to a partial differential equation (Schrodinger equation) by a finite system of coupled but 
ordinary differential equations for finitely many variables. 

For a given ordinary (unconstrained) quantum system 2 specified by a Hilbert space TC 
and a Hamiltonian H, the evolution is given as a flow on TC by the Schrodinger equation: 

-ih±Mt)) = lft\9(t)). (2.1) 

Equipped with the inner product (-|-) in TC, any operator F on TC defines a function 
F := (F) mapping G TC to the complex number (\E r |.F|\E r ). For two functions F = (F) 
and K = (K) , the inner product of TC defines the Poisson bracket of F and K: 

{F,K} = ±([F,K}), (2.2) 



1 The reader may find Reference [To| useful to serve as a quick review for the background knowledge needed. 

2 In this paper, we only consider quantum mechanical systems, not quantum field theories. 
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which again maps |\&) G TC to a complex number. For example, q := (g) and p := (p) lead 
to {q,p} = 1 from [q,p] = ih. 

Let {q\pi) be the set of fundamental operators on 7i, we define the classical variables as 

:= (<f >, Pi : = (2.3) 
and the quantum variables of the n-th order as 

Qh...i n ._ _ ^ . . . _ a^))^ = _ . . . (£<n) _ (2.4) 

associated with the fundamental operators {x*}i<j< 2 Ar := {g fe ,Pfc}i</c<Af, where the subscript 
Weyl denotes symmetric ordering of the operators as the parenthesized indices represent the 
symmetrization: A^^ = l/n\ Yip Ai P{1) ...i P{n) for all permutations P . These classical and 
quantum variables are bounded by Schwarz inequalities. A special case is the well-known 
uncertainly relation: ((q — q) 2 ) ((p — p) 2 ) > {{q — q){p — p))\f eyl + h 2 / A> h 2 /4. 

In terms of classical and quantum variables, the Schrodinger equation (12.11) govern by 
the Hamiltonian H, which is usually taken to be the Weyl ordered operator if(x l )w ey i with 
H(x l ) being the classical counterpart, the quantum evolution can be equivalently described 
as the Hamiltonian flow: 

x l = {x\H Q }, (2.5) 
Qh...i n = ^G il - in ,H Q ], (2.6) 



d n H(x^ 



generated by the quantum Hamiltonian defined as 

n=2 

The reformulation gives the equations of motion in a classical form (i.e., a system of 
ordinary differential equations, albeit for infinitely many variables) and makes it possible to 
analyze the classical limit in a direct manner with quantum effects taken into account. While 
the classical variables represent the expectation values of the wave-function, the quantum 
variables carry the additional information (e.g., spreading, squeezing, etc.) of the wave- 
function around the peak. If the wave-function is highly coherent and thus sharply peaked 
around the semiclassical trajectory, as a good approximation, we can truncate the infinite set 
of coupled equations of motion to a finite set of differential equations by ignoring the higher 
order quantum variables. This offers a systematic method to study the effective equations 
of motion order by order. (For linear systems, such as the harmonic oscillator, the infinite 
set of coupled equations decouples into many sectors, each of which contains only finitely 
many variables.) 



B. Effective equations for constrained quantum systems: Approach I 

Many quantum theories of our interest (such as general relativity and cosmology) are 
formulated as constrained quantum systems, for which, however, the effective theory devel- 
oped in [l3| and outlined in Sec. [TT1 cannot be directly applied. For a constrained system, 
the dynamical evolution is completely given by the Hamiltonian constraint equation: 

H\V) = 0, (2.8) 
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which dictates the correlations between physical variables, rather than the physical variables 
evolving with respect to a preferred time variable. Because the dynamics does not depend on 
the coordinate time explicitly, the effective theory for constrained quantum systems requires 
more considerations than that for ordinary quantum systems and its rigorous formulation is 
under development fl~8j ]. 

Nevertheless, for the constrained system that possesses an internal clock, there are ways 
to circumvent the difficulty arising from the absence of the explicit time variable, allowing us 
to analyze the dynamics in the same manner of an ordinary quantum system. To elaborate, 
let us consider a spatially homogeneous cosmological model with a scalar field <fi. By choosing 
the lapse function iV associated with the time coordinate t' (i.e. dr = Ndt' with r being 
the proper time), the Hamiltonian constraint can be rescaled by N to be of the form: 

H = Ho + ^ (2.9) 

where p^ is the conjugate momentum of (ft, satisfying {</>,£></,} = 1, and Ho is the rest part 
including the potential of (if any). The Hamiltonian constraint equation (12. 8p then yields 

Pi\V) = -ih^r\9) = ±(-2£ ) 1/2 |^> =: ±he 1 / 2 \V), (2.10) 

where the sign ± is for positive/negative frequency solutions (i.e. expanding/contracting 
solutions with respect to 0), respectively. If the positive and negative frequency sectors 
do not interfere with each other (as in the case that <fi is free of potential), we can restrict 
ourselves to one of them. Consequently, compared with the ordinary quantum system govern 
by (12. ip . the variable (j) now plays the role of the time variable t while (ft 2 0) 1//2 serves as the 
ordinary Hamiltonian. 3 That is, the quantum evolution of a constrained system has been 
reformulated as an ordinary quantum system, with the scalar field treated as the ordinary 



time variable. This is the approach adopted in 16] to study the effective equations of motion 



for the quantum cosmology in the k = FRW model. 



C. Effective equations for constrained quantum systems: Approach II 

For the quantum cosmology in the k = FRW model, the square root of remains 
polynomial of the fundamental operators and Approach I can be easily applied. Moreover, 
the Hamiltonian gives a linear system and as a consequence the equations of motions are 
exactly solvable 16| . 



However, for more complicated cases, such as the quantum cosmology in the Bianchi I 
model, which will be the main focus of the rest of this paper, 1//2 is no longer polynomial of 
the fundamental operators and because of the non-polynomiality the approximation scheme 



In general, may not be positive definite, but in the fully developed quantum theory, it turns out that 
the subspace of the negative eigenvalues of is nonphysical and thus the square root of is well-defined 



in the physical subspace on the positive spectrum of 0. This can be easily understood if we rewrite 
as f^l'S) = h 2 Q\^), which excludes the negative spectrum of since on the left-hand side is positive 
definite. More details can be found in Q for the Bicnchi I cosmology studied. 
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of Approach I is no longer under good control (as can be seen in Appendix |A|). To work 
around the technical difficulties, we innovate a second approach — Approach II. 

The idea of Approach II is to treat the constrained quantum system as if there was no 
constraint in the first place. The formulation outlined in Sec. Ill Al for ordinary quantum 
systems is then used to obtained the equations of motions with respect to the coordinate 
time t'; that is, we define Hq via (12.71) associated with H given by (12. 9p and have 

{x\H Q }, (2.11) 

{(r-".H (J \ for fc = l, 2,..- ,n, (2.12) 

as the equations of motion for the classical and quantum variables up to the n-th order 
of our interest. On top of these finite differential equations, we then weakly impose the 
Hamiltonian constraint as additional relations, by which the equations of motion abide, by 
demanding 

(y\(h 2 Q) k / 2 \V) = (M/|(±p ) fc |vI/> for k = 1, 2, ■ ■ • ,n, (2.13) 

up to the n-th order. If we only consider quantum variables up to the n-th order and neglect 
the higher order variables in (I2.12p . correspondingly, we should impose the constraint (12. 13)) 
up to the same order. The n constraints of (12. 13)) will correlate the classical and quantum 
variables of the order < n while higher order variables are consistently ignored. In the 
end, we eliminate t' in favor of <fi through the relation d<f>/dt' given in (12.111) . The resulting 
dynamics is in a form independent of the time coordinate t' and predicts only the correlation 
between physical variables and 0, which now serves as the internal time. 

Approach II, however, is not entirely equivalent to Approach I. For the WDW theory in 
the Bianchi I model, Approach II is studied in Sec. II II I and, for comparison, Approach I is 
included in Appendix |A] The comparison shows that the two approaches yield the same 
result at the 1st order but already disagree at the 2nd order. When compared with the 
results of the fully developed quantum theory, it turns out that Approach II give sensible 
results while the order-by-order approximation scheme of Approach I is messed up due to 
the non-polynomiality. 

Therefore, even though Approach II is only heuristically motivated, it gives a very good 
effective description of the quantum dynamics. If we take the viability of Approach II 
very seriously, we can even postulate Approach II as an alternative formulation for the 
fundamental theory of the constrained quantum system that possesses an internal clock. In 
order to be a viable formulation, in addition to giving the prediction very close to that of 
the standard treatment, Approach II has to be checked for self-consistency at least by two 
tests. First, for a given initial state which satisfies (12 . 13)) . after being evolved via (12. lip and 
(12. 12)) . the evolved state should satisfy (12 . 13)) again (up to the n-th order). Second, if we 
choose a different lapse function associated with a different time coordinate t", H in (12.91) is 
rescaled accordingly and so is Hq. The equations of motion given by (12. lip and (12. 12)) are 
then with the rescaled Hq and with respect to t" but they should yield the same result (up to 
the order of our interest) when related to the internal time, regardless of different choice of 
the lapse function and coordinate time. If this is case, the coordinate time is really nothing 
but an auxiliary variable. Both tests are justified in Appendix [B] but the justification is not 
completely stringent for generic cases. 

In the rest of this paper, we focus on the quantum theory of cosmology in the Bianchi 
I model, which sets a benchmark to test the reliability of Approach II. As we will see, 



dx 1 

~dV 

dG h-i k 

dt' 
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both for the WDW theory and simplified LQC, at least up to the 2nd order, the modified 
constraint (I2.13P simply correlates the constants of motion, and therefore the first test is 
trivially affirmed. Moreover, the studies of Approach II for the WDW theory and simplified 
LQC in the Bianchi I model give very sensible results compared to those obtained in the 
fully developed quantum theory. This is very instructive and suggests that Approach II, 
albeit not rigorously developed, could represent a sound quantum theory of its own even at 
the fundamental level. 

In the standard procedure to quantize the constrained system, as the constraint is imposed 
via (12.81) . the physical states that satisfy the constraint are distributional if zero is part of 
the continuous spectrum of the constraint operator H; in this case the physical states are 
not in the kinematic Hilbert space 7ikm but instead in the dual space of 7ikin- Because of 
the distributional feature, the physical Hilbert space 7i p h ys is rather difficult to construct 
and a variety of techniques have to be applied such as group averaging and refined algebraic 
quantization [l9|, 20] or the procedures performed in 0, 0] to define the physical inner 



product of Tiphys by identifying a complete set of Dirac operators. On the other hand, when 
the scheme of Approach II is performed, we approximate the quantum system by finitely 
many variables relevant to the accuracy of our measurement. Since we truncate the infinite 
degrees of freedom to finitely many variables and, accordingly, the Hamiltonian constraint is 
weakly imposed only to relate these finite degrees, the issue of the distributional property is 
completely gone and the kinematic Hilbert space remains the arena for the dynamics. The 
technical difficulties and conceptual obscurity of the standard treatment for constrained 
quantum systems are thus avoided in Approach II, whereby the evolution equations (12. lip 
and (I2.12p subject to (12.131) can be intuitively understood as the direct extension of the 
classical Hamilton's equations and the classical Hamiltonian constraint. 

Additionally, Approach II has one extra merit as measurement of the internal time is 
regarded. Both in Approach I and in the standard treatment for the constrained quantum 
system, dynamical evolution is described as the correlation between physical variables and 
the internal time 0, but, because of the aforementioned distributional property, plays a 
special role very different from other observables: In the resulting formalism, is treated 
as a pure parameter (c-number) as opposed to the other physical observables, which are 
treated as operators (q-numbers). In other words, the operator is not defined in the 
physical Hilbert space (p^ is well-defined though) and therefore there are no notions of 
measurement, uncertainty and so on for unless inferred indirectly. In this regard, the 
philosophy of "timeless" formulation for relativistic quantum mechanics as advocated in 
[2ll . 22| is not fully realized in the standard treatment for constrained quantum systems. In 



Approach II, on the other hand, is put on the same footing as all other variables. Not 
only we can define the expectation value (0), uncertainty A 2 (0) := and so on for 0, 
but the equations of motion (12.111) and (I2.12p also tell us how they evolve with time. 4 The 
variable is no distinct from other observables except that it is also used as an internal 
time to label the evolution. 5 From the timeless perspective, Approaches II seems to give 
a more satisfactory framework such that all physical observables, including the internal 



time, correspond to measurements and every measurement yields uncertainty (see also |23|). 



4 More comments about the evolution of the variables involving <fi can be found in the end of Sec. IIII CI 

5 The reason why <j), not other observables, serves as the internal time is because we take the square root in 
(|2.13p with respect to . It is explained in Appendix [B] that the free massless scalar is special for being 
the internal time as the first test is concerned. 
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Compared with (12.81) for the standard treatment and f 1 2 . 1 j) for Approach I, the notable 
difference of Approach II is that the Hamiltonian constraint is only weakly imposed as (I2.13P 
consistently ignores quantum fluctuations of higher orders. It seems to be the ignorance that 
retains the notions of measurement and uncertainty for 0. This makes us suspect that the 
quantum Hamiltonian constraint (12.81) may turn out too restrictive and nonphysical, as it 
relates quantum fluctuations to every level of depth, which could be irrelevant after all as 
far as the physical measurements with intrinsic uncertainties are concerned. 

Although some implications of Approach II remain speculative, the study of it can still 
teaches us valuable lessons for solving constrained quantum systems. In the following, we 
will apply Approach II to study the WDW theory and simplified LQC in the Bianchi I model 
as two sound examples and see explicitly how the ideas of Approach II materialize. 



III. EFFECTIVE EQUATIONS OF MOTION FOR THE WDW THEORY 

With the formulation of Approach II at hand, we are ready to study the effective equa- 
tions of motion for the WDW theory in the Bianchi I model. The results, except those 
regarding measurement of <ft, will be shown to agree very well with those obtained from the 
fully developed quantum theory and thus the viability of Approach II is affirmed. On the 
other hand, Approach I is presented in Appendix |A] for the WDW theory and shown to be 
incompetent for the Bianchi I model. 



A. WDW theory in the Bianchi I model 

The phase space of the Bianchi I model is given by the diagonal triad variables pi and 
connection variables c 1 for I = 1,2, 3. The diagonal variables satisfy the canonical relation: 

{c I , Pj } = 87iG 1 6 I J . (3.1) 

In the presence of a free massless scalar field <fr(x,t) = <p(t), (which is independent of the 
spatial coordinates with homogeneity assumed), the classical dynamics is govern by the 
Hamiltonian constraint: 

H = ifgrav + 

1 £> 2 

= - „ 2 - = (c 2 c 3 p 2 p 3 + cVpiPa + c l c 2 Pl p 2 ) + —A , (3.2) 

where p^ is the conjugate momentum of <j) and has the canonical relation: 

= 1. (3-3) 

We can rescale the Hamiltonian by choosing the lapse function N = y/Pip 2 P3 associated 
with the new time variable t' via dt' = (pip 2 P3y 1/2 dT. The new (rescaled) Hamiltonian 
constraint reads as 

H = ~1T7^ (<?(?Pm + c'cWs + cVp^) + ^. (3.4) 
87rG7 z 2 
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In the Wheeler-DeWitt (WDW) theory, the standard Schrodinger quantization is 
adopted. The kinematic Hilbert space is T~C^ W = L 2 (R 3 , d 3 p) ®L 2 (IR, d(f)). The variables c 1 , 
Pi, 4>, p^ are lifted to the operators c 1 , pi, 4> and p^. In particular, &i and p^ are represented 
in 7iSrP w as differential operators: c 1 — > ih&TtG^-J^, p^ — > —ih-^. We also define the 
hermitian operators: 

^ / = ^(c 7 P/+P/c 7 ) (3.5) 
and correspondingly define the Hamiltonian operator as 

# = (^ 3 + ^ + ^K 2 ) + f =: -\ {n 2 6 - p 2 } , (3.6) 

by fixing the ordering such that H is hermitian. 6 



B. Classical and quantum variables 



In order to apply Approach II to derive the effective equations of motion, we define 
classical variables as 7 



Pi 



(Pi), 
{P<f>)> 



K 1 := (K 1 ), 
:= (4>), 



which, by (12.21) . satisfy the Poisson brackets: 

{p h K J } = -torGySfa, {0, p*} = 1, { P i, P j} = {K 1 , K J } 
The associated quantum variables of the 2nd order are also defined: 

G (n=2) . 



QUn=2) 



Gn 





= ((Pi ' 


- (Pi))(pj ~ (Pj))), 


G IJ 


= ((K 1 


- (&)){&-{&))), 




= ((K 1 


- {K'))(pj - 


P^Pcf, 


= ((pV 


- (P<t>))(P<p- (p*))), 




= ((#- 




G^ 

p<i> 


= ((#- 


(4>))(P<t> - (P»)) W eyl- 



(3.7) 
(3.8) 



(3.9) 



(3.10) 



(3.11) 



We have {(j),G^ n 2 )} = {p^,G^^ n 2 )} = 0, but on the other hand, as opposed to the case 
for canonical variables (p/,c J ), the commutator \pj,K J ] is not a constant and consequently 



6 The WDW theory constructed here is very similar to that in [8(. In particular, 8«6 (9 defined in [8() 
with very tiny difference due to the different choices of ordering. is not positive definite but the square 
root of 6 is well-defined in the physical subspace. Also see Footnote El 

7 For technical convenience, instead of the canonical variables (pi,c J ), we choose (pi,K J ) as the fun- 
damental variables. For reference, the Poisson brackets based on the canonical variables are listed in 
Appendix [Cl 
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the quantum variables do not commute with the classical variables but, by (12.21) . yield the 
relations: 

{pi, G JK } = 0, { PI , G JK } = -8ttG 7 + 5fG J K ) , 

{ P i,G J K } = -SttG^Gjk, (3.12) 

and 

{K 1 , G JK } = -8ttG 7 {S r j + S T K ) G JK , {K 1 , G JK } = 0, 

{K 1 , G J K ] = -MG^Gi. (3.13) 

We also have 

{G**, } = 2G**, {G**, G Pm } = 4G* , {G* , G HH } = 2G W , (3.14) 

and 

{G J; Gf } = (8vrG 7 ) {SiCfij - 5«Gi K ) , {G IJ , G KL } = 4(8vrG 7 )^G^, 

{G IJ , G K L } = 2(8^7)4^^, {Gu, Gl } = -2(87tG 7 )c^G j)L k, 

{G u , G KL } = {G /J , G KL } = 0, (3.15) 

where the 3rd order quantum variables (G*™^) are defined in the obvious manner (e.g., 

Gj K '■= ({K 1 - (k I ))(p J - (pj))(pK - (pK))) Wey i)- The Schwartz inequality leads to the 
relations: 

G*+G p ^>G$; + ^ (3.16) 

G n G JJ > G IJ \ (3.17) 

GnGjj > Gu 2 , (3.18) 

G"G„ > G? + \\(IK;PA)\ 2 = Gj 2 + ^Is'pl (3.19) 



and 



C. Effective equations of motion 

Corresponding to (13.61) . the quantum Hamiltonian defined in (12. 7p reads as 
Hq = -g^^ (K 2 K 3 + K X K 3 + K X K 2 + G 23 + G 13 + G 12 ) + ^ + ^p. (3.20) 

According to (12.111) and the Poisson brackets listed in Sec. IIII B\ the equations of motion 
for the classical variables are given by 

dK 1 



dt' 
dpi 

dt' 
d4 

dt' 

dpcf, 

dt' 



{K 1 , H Q } = K 1 := 87rG 7 fr/C 7 are constant, (3.21) 

{ Pu H Q } = r l {(K 2 + K 3 ) Pl + G 2 + G 3 } , (3.22) 
{<t>,H Q }=p h (3.23) 



{0, Hq} = =>- Pcj, '■= ^VovrG/C^, is constant, (3.24) 
11 



where we define the dimensionless constants KJ and K,^ for convenience. We also have the 
equations for dp2/dt' and dp^jdi! in the cyclic rearrangement of (I3.22p . 8 Similarly, (I2.12p 
gives the equations of motion for the quantum variables: 

dG IJ 

= {G IJ ,H Q } = => G IJ are constant, (3.25) 



dV 
dG\2 

dt' 
dd 



{G 12 , H Q ] = 7~ 1 G 12 (K 1 + K 2 + 2K 3 ) + 7 - 1 0(G< n=3 >), (3.26) 



w l - = {G u , H Q } = T - 1 Gn {2K 2 + 2K 3 ) + ^(G^), (3.27) 
^ = {G{, H Q } = r l G[ (K 2 + K 3 ) + j^O(G^% (3.28) 



and 



dG 



p<t,p<t, 



{ G pm, h q}=° G 9 . p . is constant, (3.29) 



dt' 

—^r = \ G l^ h q} = g pop<p' ( 3 - 30 ) 

dG^ 

= {G»,H Q } = 2G^. (3.31) 

Note that f!3.23j) and (I3.24p imply that is a monotonic function of t' and thus can 
be used as an internal time. If the wave-function is sharply peaked, (9(G <n=3 - ) ) terms are 
negligible, and, in terms of the internal time 4>, the general solutions to the above equations 
of motion are given by 

G /J (0) = g IJ , (3.32) 



GM) ~ 9ue 



(3.33) 



/ 87rG( ^cM )( ^ o) 



Gn(0) « Pne V ^ (3.34) 

r/ , T Vs^gI )u-<p ) 
G{(4>) ~ ffi e V 5 ^ (3.35) 



and 



/,(-) = (p liO + ^^-g -(0-0o)l> v f ; • (3.3G) 



G p^p^{4>) — 9v$pv> (3.37) 



M = < + H==7^ - (3-38) 



hV&xG^ 
2q* 



gt"t> + - (0 - O ) + - I^Z. , ( o - oo ) ^ . ( 3.39 i 



The obvious cyclic repetition will not be mentioned again hereafter. 
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where pi t o, g IJ , gu, g^, etc. are all constants, which satisfy the inequalities 

g n 9 JJ >g Ij2 , 9ii9jj>9ij 2 , for V I, J, (3.40) 

g H gjj>gj 2 , for I ^ J, (3.41) 

11 ^ 12 . (87TG/17) 2 / £ 2 + # 3 \ 

5 9n>9i + - A Pi,o + —^==—(0-00) , (3.42) 



and 



according to (l3~T6l)-(l3~T9j) . 

Later in Sec. IIIIDl we will show that the Hamiltonian constraint gives the rela- 
tion (13.471) for the dimensionless constants KJ and JC^, which agrees with the classical 
counterpart with small corrections. Consequently, compared with the classical solution: 
Pi(0) = pi j0 exp(vo7rG(0 — O )(/C 2 + /C 3 )//^), (13.361) shows that the expectation value of 
each pi closely follows the classical trajectory but will receive back- reaction if the off-diagonal 
squeezing is nonzero (i.e. gf,gf 7^ etc). 9 Furthermore, in the case without off-diagonal 
squeezing (i.e., gf,gf = etc), according to f!3.34|) and (13.361) . the relative spread for each 
"pi remains constant: i.e., Api/pi := y/Gji/pi = y/gli/Pi^o- 

The fully developed quantum theory of the Bianchi I WDW cosmology has been studied 
in Section VI and Appendix A of [8]. The semiclassical state constructed in (A.l) of [s| 
is with constant spreads in Ki (i.e. G n = (87i , yh) 2 a'j), zero off-diagonal correlations (i.e. 
Gu = for / 7^ J) and zero squeezing (i.e. Gj = 0). The evolution of this state is given by 
Equations (6.19), (6.21) and (6.22) in [8], which indicate that the peak of the wave-function 
follows the classical trajectory with constant relative spreads. That said, the behaviors of 
f l3.32p - fl3.36p conform to the fully developed quantum theory at least for the case studied in 
[8|. The results here also instructs us to construct a new type of coherent states with nonzero 
squeezing and off-diagonal correlations by generalizing (A.l) in [8[ to see if the evolution in 
the fully developed quantum theory receives the back- reaction as suggested in f)3.36p . We 
expect the answer to be affirmative. By contrast, as shown in Appendix |Al the scheme of 
Approach I does not give the same results as those in ^ and is thus inept for the Bianchi I 
model. 

Meanwhile, as mentioned in Sec. Ill C\ the notions for the measurements involving are 
well-posed in Approach II and furthermore how they evolve can be answered. According to 
(I3.37p - fl3.39p . the spread of p^ remains constant, but the squeezing and the spread of grow 
with 0. This is a well-known fact as for the quantum mechanics of a free particle, which 
has a strictly growing spread. For unbounded systems, we cannot expect to have a valid 
semiclassical approximation for all times, but the semiclassical treatment is still reasonable 



for limited amounts of time (see Example 3 in [13j for more details). This implies that 
treating as a classical clock is only valid for a limited period of time. This problem 
can be mended by adding a mass term to 0, which then behaves as a harmonic oscillator 
instead of a free particle. In the presence of the mass term, as for a harmonic oscillator, the 



9 In this paper, squeezing refers to the correlation between pi and K (i.e. Gf) or between cf> and (i.e. 
Gp ). The phrase "off-diagonal" refers to Gu, G IJ and Gj with I ^ J while "diagonal" to those with 
I = J. 
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wave-function always admits bounded spread and thus the semiclassicality can be sustained 
forever. Even though is no longer a monotonic function of t', it can still be used locally 
as a classical clock for all times. It seems to suggests that only the field with oscillatory 
behavior of harmonic type can be viewed as a classical clock. This observation might shed 
light on the problem of time (see 24j for an overview) by offering insight into why in reality 
the passage of time is always measured in terms of oscillatory signals. 



D. Hamiltonian constraint 

In addition to solving the effective equations of motion, we still have to impose the 
Hamiltonian constraint (I2.13p . which will give further relations to relate the constants of 
motion used to parameterize the semiclassical solutions. Since we have solved the equations 
of motion up to the order of G^ n=2 ^ and G^^ n=2 \ accordingly, we only consider k = 1,2 in 
(I2.13p . (The Hamiltonian constraint for k > 3 will further constrain the higher order terms 
Q{<t>,n>2,) Q(n>3) anc } Q^gj. parameter constants.) 

For k = 2, fl2TT3D yields 

(Pi) = ^{(K^ + iK'K^ + iK'K 2 )} 

= {K 2 K 3 + K X K 3 + K X K 2 + G 23 + G 13 + G 12 } 

= 16nGh {1C 2 1C 3 + K 1 1C 3 + IC'K 2 } + -^- 2 {g 23 + g 13 + g 12 } . (3.44) 

For k= 1, by defining O^ 1 ) := (K 2 K 3 + K 1 K 3 + K 1 K 2 ) and y/6(K J ) := (K 2 K 3 + K 1 K 3 + 
K 1 ^) 1 / 2 , (12331 yields 

= 7th^ (K ' + {K ' - K ' ))} 

d n Ve 



if 00 l 

| v^') + g - dKh _ _ _ dKIn 



G 



hh-In 



1 J ^ - + **>' G n _ + **)' g 22 - (f ' +£) 2 G 33 ( 3. 45) 



v/47tG7 2 1 8^(^) 3 / 2 8#(if 7 )3/ 2 8#(if 7 )3/ 2 

, djK 1 ) - k 12 r23 ejK 1 ) - k 22 ejK 1 ) - k 32 ^- n n<r^y 

+ 49(K*)W ° + 4 5(^)3/2 6 + 4(9(^)3/2 6 " > 



and consequently 
1 



„ > W ) _2W tfl _... + ? m^ G . + ... i ( , 46) 
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which can be rewritten as 



(/C 2 /C 3 + K}K? + KMC 2 ) - K> , : , 



2(/C 2 /C 3 + /OK? + /G/C 2 ) 



9 M + ■■■>, (3.47) 



closely in agreement with the classical counterpart /C| = 2 (/C 2 /C 3 + /C^/C 3 + K}K?) with 
small corrections arising from the quantum fluctuations. 

Furthermore, (13.441) and (I3.45j) together give the spread of p$. 

(Ap ) 2 := (pj) - (ptf,) 2 = G nV4> = g nn 



(47TG7 2 )- 1 f 
(/C 2 /C 3 + /C 1 /C 3 + /C 1 /C 2 )\4 



1 [(k 2 + ^y 1 + (k 1 + k 3 )V 2 + (k 1 + k 2 )V 3 ] 



i 

+ 2 



K'V 3 + K 2 V 3 + K 3 V 2 ] | + < - iwG '^ 1 ( 9 » + 9 i3 + 9 u) (3 48) 

In the case without off-diagonal spreads in i^ 7 (i.e. g IJ = for I ^ J), this leads to 



« pi {(^ 2 + K 3 )V 2 + (K 1 + K?fa\ + (K 1 + K 2 )V 3 2 } , (3.49) 

where the dimensionless constants a 2 := G I ^J(87rG , yh) 2 are the spreads for K 1 /(8irG , yh). 
This result concurs with Equation (6.24) of |8| for the WDW semiclassical state except the 
mismatch of an overall factor 2 due to the approximation method used in jsf to evaluate 
Ap^. The agreement further strengthens the viability of Approach II. 

IV. EFFECTIVE EQUATIONS OF MOTION FOR SIMPLIFIED LQC 

The viability of Approach II has been tested for the WDW theory. In this section, we 
will investigate the effective equations of motion for the simplified LQC in the Bianchi I 
model along the line of Approach II. The results are expected to be in accordance with 
those obtained in the semiclassical approach [HI 12] amended with back- reaction resulting 
from quantum fluctuations. 



A. Simplified LQC in the Bianchi I model 

The detailed construction for the Hamiltonian operator corresponding to (13. 2p in the 
fully developed LQC can be found in The LQC quantization is implemented with two 
main sources of quantum corrections: First, the connection variables c 1 do not exist and 

— 1/2 

should be replace by holonomies; second, the inverse triads p l , upon quantization, receive 
quantum corrections via the so-called Thiemann's trick. It is realized that the corrections 
on the inverse triads are less important when the wave-function evolves in the semiclassical 
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realm. To simplify the quantum theory, the Thiemann's trick is simply ignored; i.e., the 
Hamiltonian constraint (13.21) is trivially rescaled to (13.41) before quantization. 10 

Therefore, we start with the rescaled Hamiltonian (13. 4p and take the prescription to 
replace c 1 by sin(/2/c 7 ) / \xj by introducing discreteness variables fij. The rescaled Hamiltonian 
is modified to 

- 1 Jsin(^2C 2 )sin(/i3C 3 ) \ t p l fA i\ 

H = o r< <? i — P2P3 + cyclic terms > + — . (4.1) 

There is a variety of possibilities to implement the fx discreteness and two well-motivated 
constructions (referred to as the "/I scheme" and il p! scheme'') arc discussed in great detail 
in [12|. The fl' scheme has a better scaling property that the quantum dynamics based on 
it is invariant under different choice of the fiducial cell. Unfortunately, the quantum theory 
of the fl' scheme is much difficult and has yet to be constructed. In this paper, we study 
the /z-scheme quantization described in [s| with the aforementioned simplification. In the 
/2-scheme, fij are prescribed as 






Pi = \—> t*2=\—, /i 3 =W— , (4.2) 



where A is the area gap in the full theory of LQG and A = 2 (^^pi) ^ or ^ ne standard 
choice (but other choices are also possible) with £p\ = \J Gh being the Planck length. 

The kinematic Hilbert space for LQC is given by ftj£ c = L 2 (R| ohr , d 3 p Bohv ) <g> L 2 (R, dcj>), 
where <fi is in the ordinary Schrodinger representation while p is in the "polymer represen- 
tation" , which reflects the very feature of LQC that the connection operators c 1 cannot be 
defined. Our task now is to promote (14. ip to a well-defined quantum operator in Wjjvjj . 
Following the strategy used in [8], we define the dimensionless affine variables vj to satisfy 

^^•wrh (4 - 3) 



which yields 



3/2 
Pi Pi 



pi ~ A i /2 " 6tt 7 4i^. (4.4) 

In L 2 (R| ohr , d 3 pBohr) = -^ 2 (^Bohr; ^^Bohr), we can then quantize exp(±i/2/c / /2) as the trans- 
lation operators 

e ±2Mic ._ e +a n ^(#) — ^( Vl rp 1,^3,^3) (4.5) 

and so on in the ^-representation. 11 Here, the wave function ip{v) is related to ip(p) via 

/ A 3 \ 1/4 

\ijj(v)\ 2 d 3 v = \^(p)\ 2 d 3 p => ip(v) = (4irp£ 2 ) 3/2 [ </>(£). (4.6) 

\P1P2P3J 



10 This the simplified LQC studied in [J,|5| for the isotropic case, where the simplification leads to the exact 
solvability. 

11 If the Hilbert space was in the ordinary Schrodinger representation, c 1 would be well-defined and rep- 
resented as c 1 — > ift87rG7g|j in the ^-representation. In the ordinary quantum mechanics, p — > — jr-^ 
and e la P/ h ip(x) — ij)(x + a), but we cannot simply take e l ^ x ^/ n — > e^ x ^ because it is not unitary in 
general. Here, by contrast, since p is in the polymer representation, the same problem does not occur. 
The operators defined in (|4.5|) are unitary in L 2 (R| ohr , d 3 t>Bohr)- (Also see Footnote 10 in Q). 
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In L 2 (M Bohr , (i^Bohr), we then define the unitary operators 

j 1 ■= e~ 2 A ~ e +iilicl or JV(v) = ip{vi - 2, v 2 , v 3 ), (4.7) 
J 7 t := e + 2 ^7 ~ e"^ icl or J 1 ^^) = + i*). (4.8) 



Note that J 1 and J 7 "!" correspond to their classical counterparts e +iflicl and e *^ icl ; respec- 
tively, and we have 

[v z , J J ] = 25j J J , [uj, J Jt ] = -25jj J \ J 7 J Jt = J Jt J 7 = (4.9) 
For later use, We also define the hermitian operators 

K 1 : = 67T74, — fej 7 - J 7t ^l , (4.10) 
2% L J 

£ 7 := 67r 7 £ P1 i [^J 7 + J 7t f) 7 ] , (4.11) 
which correspond to the classical counterparts: 

K 1 ~ si ^ 7C \ — P/ c 7 , (4.12) 

/X/ /Ijc'->0 

^ „ C0 <fii cI ) pi = 67r7 4 lC os(/2 lC 7 )^ — 6^74^/ = ^ /2 A- 1/2 . (4.13) 

/i/ /tire*— >0 

By ( 1431) . we have 

i^ J ] = -2i5jSj, [v T , Sj] = 2t5 I jK I } [K 1 , Sj] = -21(6^^% (1 - «j) , (4-14) 
and obviously 

[^,{; J ] = [K 7 ,^ J ] = [^,^ J ]=0. (4.15) 

Now, in the kinematic Hilbert space to be TC]^ = L 2 (IR Bohr , <i 3 i>Bohr) ®L 2 (R, d<f)), let the 
Hamiltonian given by (14.11) promoted to the hermitian Hamiltonian operator: 

H = --4— (k 2 K 3 + K'K 3 + K'K 2 ) + % (4.16) 

where — > —ifr-jk is the same as that of the WDW theory studied in Sec. IIHI The quantum 
Hamiltonian given by ( 14.161) is exactly the same as that of the simplified LQC studied in 
We will take (14.161) as our fundamental quantum theory and study its effective equations of 
motion, which are to be compared with the numerical results obtained in [d]. 12 



12 Both in the detailed construction of LQC [8| and in the simplified LQC [9|, to extract the physical 
information, we have to further restrict ourselves to one of the super-selected sectors 7i phys as our physical 
Hilbert space and define the physical inner product ("l^p^ys on it- In Approach II, as commented in 
Sec. Ill CI we do not need to construct the physical Hilbert space Tt p h ys but still all the information about 
dynamics can be extracted solely from the kinematic Hilbert space 7ikin endowed with the original inner 
product (-|-). 
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B. Classical and quantum variables 



Before we tackle the effective equations of motion for the simplified LQC, we have to 
define the classical and quantum variables and the Poisson brackets of them. By choosing 
(vj, K J ) or (Si, K j ) as the fundamental variables for the gravitational part, the classical 
variables are defined as 



vi:=(vi), K'-^iK 1 ), S r .= (Si) 
in addition to (13.81) . By (I4.14p and f)4.15p . (I2.2p then gives the Poisson brackets: 

2 



and 



{v I ,K J } = ~6jS I , {v I ,Sj} = ^5 I J K I , {K'.S.,} jiS- 



{vi,v J } = {K I ,K J } = {Si,Sj}=0. 



-vj) 



(4.17) 

(4.18) 
(4.19) 



The associated quantum variables of the 2nd order are also defined: 



£r(n=2) . 



Gu 
G IJ 

and Gpji 
Gj> 



((Si-(Si))(Sj-(Sj))), 
((k I -(R I ))(R J -(R J ))), 

((k 1 - (k 1 ))^ - (Sj))) Weyl , 

(67ry4i) a <(*/-<*/»(*J-<*J»>, 
(6n 1 el^(S I - (^(vj - (vj))) Weyl , 

(67T 7 4 1 )((^-(^))(^-(^))) W e y l- 



(4.20) 



(4.21) 



Here, the primed lower indices are used for G^ 11 2 -* whenever vi is referred to; the unprimed 

lower indices are reserved for Si. 

Note that vi, K 1 and Si are not independent of one another but related via 



a, 2 -2 

K 1 +Si 



(6irjl 



ll)2 [(vi + 2f+v*]. 



(4.22) 



We suppose that the universe does not go down to the deep Planckian regime, as the big 
bounces of vi are expected to takes place at a much larger scale. 13 Hence, we assume Vi 1 
and then have 

K l2 + Si 2 ss (6n-f£ 2 Pl ) 2 vj - G H - G n + G vv . (4.23) 
For the same reason, up to the error due to quantum ordering, we can approximate 



Gu pa cosmic 1 ) cos(njc J )Grji + 
Gu, pa co S (fiic I )G r j, + 0(tl, 1 ), 
Gj ^ cos(/i JC J )G^ +0(40 ■ 



(4.24) 



13 As will be seen in (|4.52j) . this assumption can always be achieved by increasing (p<j>). This also justifies 
why the quantum corrections on the inverse triad due to Thiemann's trick can be ignored in the simplified 
LQC. 
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Equation (I4.23[) gives us 



iY /2 + G J/ 



(67T7£ 



2 ^2 ? ,2 
PlJ U 7 



+ 



G 



J7 



(6^7^ 



p,)M 



« ±(67T7£p 1 )f/ 

where in the second line we have used 

Gii 11 — On ~ [1 — cos 2 (/i/c / )]G/// 

(D(G( n = 2 ^ - 



( _ if /2 + G /7 iY /2 GW 



(6^74!)^/ (6tt74i)M 



1/2 



1/2 



(4.25) 



S 2 



(6^7^ 



2 ^2 



PI 



7 J 



Gri 



IT 



K 12 



(6ttj£ 



2 ^2 ? ,2 
Pi) U I 



Gr 



(4.26) 



to trade Gu for Grp. Equations (14.251) can be understood from (14.131) : While Si ~ 
(67T7£ 2 5 1 )f/ in the classical regime (fiic 1 — > and cos(/i/c / ) ~ 1), in the quantum regime 
where the LQC corrections become significant, cos(/2/c 7 ) eventually crosses zero and thus 
Si flips signs (±) when vj approaches a certain critical value. Similarly, (I4.24p leads to 



G 



i j 



± 



K j2 + G u K l2 Gvv 

1 - — I^TT^T + 



(6^7^!)^ (6^7^)4^ 



1^ 



X ± 



K j2 + G H K j2 Gjiji 

1 " ^TT^ + 



(6tt 7 4i)M (6vr 7 4i) 4 ^ 



1^ 



G 



u> 



± 



± 



K l2 + G H 



+ 



PlJ y / 



(G^^j) 2 !'! ' (67cy£ 
K j2 + G JJ K j2 Gjiji 



1/2 



G>, 



1/2 



Gj/. 



(4.27) 

(4.28) 
(4.29) 



As in Sec. HUB} the quantum variables do not commute with the classical variables, which 
are not chosen to be canonical variables. By (14.141) . (14.151) and (12. 2p . we list the Poisson 
brackets needed for the later use: 



K G JK } = ~ (5jG J K + 5fG«) , {vi, G JK } = -l (5/Gf + 5fG J K ) 



{vi,G J K } 



h 

2 

h 



{8 J iG JK - 5fG JK ) 



(4.30) 



and 



{K 1 , G JK } = 12vrG 7 (<^GW + ^Gjk) G jk , 
{K 1 , Gji K i} = YI-kG^ (5jG JK i + S^Gjik) G jk , 



{K I ,G J K } = l2nG 1 5 I K G J K „ 
{K I ,G J K ,} = l2nG 1 5 I K G J K , 



{K\G JK } = 



(4.31) 
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Similar to (j3.15p . we also have 

{G {n=2) ,G (n=2) } = 87iG-f G {n=3) . 
Finally, the Schwartz inequality gives the relations: 



QIIQJJ > & J2 



GuGjj > Gu , 

I 



G n Gjj > g 1 ; 



GppGjiji > Gpji 

2 



([K I ,Sj])' = G I J 2 + (6^)^(1 -vi) 2 , 

([K^vj]}' = G I / + (67r 1 i 2 Pl ) 2 6 I J S]. 



4 

(6tt7^ 2 - ■ ^ 



(4.32) 

(4.33) 
(4.34) 

(4.35) 
(4.36) 



Also note that all the equations regarding and p^ are exactly the same as those appeared 
in Sec. IfflBl 



C. Effective equations of motion and Hamiltonian constraint 

Corresponding to (I4.16p . the quantum Hamiltonian defined in (12 .7p reads as 
1 



Q 



(R 2 K 3 + K X K 3 + K X K 2 + G 23 + G 13 + G 12 ) + ^ + 



G 



8-6'- 2 v ' 2 2 



(4.37) 



According to (12.111) and the Poisson brackets listed in Sec. IIVB| the equations of motion 
for the classical variables are given by 



dK 1 

dt' 
dvi 



{K I ,H Q } = 



dt' {vi ' HQ} = ^Gh^ 



K := %TiG^hlC are constant, 
{{K 2 + K 3 )S 1 + Gj + G 3 }. 



(4.38) 
(4.39) 



Meanwhile, d(p/dt' and dp^/dt' are unchanged, as given by (13.231) and (13.241) . respectively. 
Similarly, by (12.121) . the quantum variables satisfy the following equations of motion: 



dG IJ 

dt' 
dGi'2' 
dt' 

dG\i\i 
dt 

dGy 



{G IJ ,H Q } = 



= =>- G IJ = g IJ are constant, 


(4.40) 


= IT 1 {Civ (K 2 + K 3 ) + G 2V (R 1 + K 3 ) } 




+ 7 - 1 C(G (n=3) ), 


(4.41) 


= -j- l G lv (2K 2 + 2K 3 ) + r l O{& n = 3 ^), 


(4.42) 


= §7 _1 <3{ (K 2 + K 3 ) +1 - l O{& n =V), 


(4.43) 



in addition to (I3T29D . (13301) and fEOB . 

The right-hand side of the above equations of motion involves Sj and G^ n=2 ^ with un- 
primed lower indices, although we consider only the time derivative of vj and G^ n=2 ^ with 
primed lower indices on the left-hand side. Using (I4.25P and (14 . 2 7 j) (14 . 2 9 [) to eliminate Sj 
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and G^ n ~ 2 ^ with unprimed lower indices in favor of Vj and G^ n ~ 2 ^ with only primed lower 
indices, we get 



Pr 



K 1 
G IJ 



SnG'yhK, 1 are constant, 



JJ 



g are constant, 



(4.44) 
(4.45) 



±- 7 ~V 



K 1 +G 



11 



K 1 Gyy 



X 



G 2 / -\- Gy 



1/2 



I (K 2 + K 3 ) + 



6-KjlpyVl 



(4.46) 



and 



P4>- 



dGi'2' 



-7 1 Gi>2>- 



± 



± 



K l2 + G 11 



+ 



K l2 G vv 

2 U„4 



1 1/2 



(if 2 + i? 3 ) 



^22 + £22 R 22 G 2 , 2 , 



1 1/2 



+ 



(GTijil^v 2 (6vr 7 4i)M 



(K 1 + K 3 ) 



dG\i\i 



dGy 



±-7-'G F1 , 



K l2 + G u 



+ 



K l2 G vv 
{<o^ 2 x ) 2 vl ' (6^t7|0M 

K l2 + G 11 K l2 G vv 
{<o^i 2 x ) 2 v 2 + (6^t40V 



1 1/2 



(2K 2 + 2K 3 ) 



1/2 



(ir 2 + K 3 ) 



(4.47) 
(4.48) 
(4.49) 



where we have used (13.231) and treat <ft as the internal time. 

When vj approaches a certain critical value, the factor ±[— ■ ■ -] 1 / 2 flips signs and conse- 
quently, by (14.461) . vi undergoes the bouncing behavior at the critical point, affirming the 
bouncing scenario predicted in [§, [Io|, H . 

If the off-diagonal squeezing is negligible (i.e. G I J ,/(67i'jip l K I vj) <C 1 for / 7^ J), the 
second term in the curly bracket of (14.461) can be dropped. Equations (I4.46P and (14.481) then 
yield 



dG r , 



dvi vi 
where the dimensionless quantity 



IV 1 



constant, 



Tj 



'•I 



A 2 t> 7 G F p 



(4.50) 



(4.51) 



is the relative spread of vj. In the case without off-diagonal squeezing, the critical value of 
vi is given by the constants K 1 , G n and 77^ as 



"^crit 



K 12 (1 + < + V 2 RI ) 



where the dimensionless constant 



vh - 



A 2 K I 



G 



11 



K 



1 2 



K 



y2 



(4.52) 



(4.53) 
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is the relative spread of K 1 . In terms of pi via (14.41) . (I4.46P can be recast as 

pi d<p 

where we define the directional density for the direction / as 





\l 61 1 


V K>4> J 


Ql, crit . 



(4.54) 



and its critical value is given by 

Ql crit 



ei := j| (4.55) 

(^) 3 (1 + *+^) _1 Pp. (4-56) 
with 

Ppi := (8ttG'7 2 A)- 1 (4.57) 

being the Planckian density. It tells that the big bounces occur up to three times, once in 
each direction of pi, whenever each of the directional density approaches its critical value 
and thus flip s the sign of (14.541) . This closely agrees with the semiclassical description in 
13 and Q except that the critical value Qi, ciit is now slightly modified by the overall 
factor (1 + rjl + ffp.j) involving the relative spreads and giving rise to the back- reaction. 
The back-reaction due to uncertainties in K 1 and vj is expected and desired, because one 
would otherwise be puzzled why the critical value £?j, C rit is exactly determined by the ratio 
of K,i and /C^, which are the parameters describing the classical behavior and should not 
completely dictate the quantum physics of the bouncing scenario. 

In the case with appreciable off-diagonal squeezing (Gj, -yt. 67i'j£p l viK J for I ^ J), the 
terms involving Gj, in (14.461) will further deviate the trajectory off / and the relative spreads 
rj^ will no longer stay constant. Furthermore, the evolution of vj in different diagonal 
directions couples to one another through the off-diagonal squeezing. Also note that the 
second term inside the curly bracket in (14.461) is oc vj 1 and therefore this back- reaction due 
to off-diagonal squeezing is most significant in the vicinity of the big bounce and negligible 
in the classical regime. 

Finally, the Hamiltonian constraint is exactly the same as that studied in Sec. lIIIDl except 
that all K 1 , K 1 and G IJ should be replaced by K 1 , K and G IJ . The validity of treating <p 
as the internal time follows the same argument given in Sec. IIII CI 



V. SYMMETRY REDUCTION TO THE ISOTROPIC MODEL 

One of the most challenging problems in LQC is to find a systematic procedure to derive 
LQC as a symmetry-reduced theory from the full theory of LQG. To shed some light on 
this issue, we can ask a similar but much easier question: "How does isotropy emerge from 
the anisotropic Bianchi I description?" With the effective theory at hand, we try to answer 
this question and make sense of the isotropy reduction from the Bianchi I model. We will 
discuss two ways of reduction: the formal and physical prescriptions. 

In the framework of effective equations of motion, similar to the Bianchi I model, the 
k = FRW model has the classical variables: 

v:=(v), K:=(k), S := (S), (5.1) 
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and the associated quantum variables are given by 

Q{n=2) . 



G(2) 

G^ 



G 



(i) 
(i) 



and G(2") 
G 



(HO 
G (i0 - 



<(£-<S>)(S-<S>)>, 
<(^-<K>)(K-(K>)>, 
((l-(^))(|-(|))) Weyl , 

(67ry4i) 2 <(« -<«»(« -<«»>, 
(67r 7 4i)(0-0»(^-^»)we yl > 

(67T 7 4l)((^ -<«»> Weyl- 



(5.2) 



(5.3) 



To get the effective equations of motion for the isotropic case from those for the Bianchi 
I model, we can formally impose the condition on the Bianchi I variables to demand that 
the classical and quantum variables with different indices are all identical; i.e., 

^1=^2=^3= V, (i.e. Pl=P2=P3= p), 

R 1 = K 2 = K 3 = K, Si = S 2 = S 3 = S, 
Gn = G22 = <5 33 = G12 = G 13 = G23 = (5(2) and so on for other G (n=2) , (5.4) 

which is trivially compatible with (14 . 3 8 [) ( 14 . 49 }) . With this condition imposed, the equations 
of motion reduce to the isotropic counterparts. In particular, (14. 46ft becomes 



dV ~ _I_Q -1 

j— « ±37 v 



n l/2 



l^ 2 



K 



where 



K = RnG^hK, and p^ = hv&nGlCq 
and the Hamiltonian constraint (13.471) now reads as 

K-l = 6JC 2 . 

Equations (15.51) can be recast as [cf. (I4.54p - (l4.56p ] 



(5.5) 

(5.6) 
(5.7) 



1 dp /16ttG 
dpd<f) V 3 

where ps is the matter density of 0: 





1/2 r 






Pcrit _ 


r 



and its critical value is given by 

Pcrit = 3 (l + rf v + rig) 1 Ppi < 3p Pi . 



(5.8) 



(5.9) 



(5.10) 



This perfectly concurs with the exact solution obtained in 0] , which indicates that the matter 
density has an absolute upper bound and the more semiclassical the state is (correspondingly, 
r]l and rj^ are smaller), the closer the upper bound can be reached at the big bounce. 
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Although this procedure formally gives rise to the isotropic theory, in the context of the 
Bianchi I model, however, the solution satisfying condition (15.41) is not "isotropic" at all. 
In fact, it is the most anisotropic one in the sense that three different directions are now 
completely correlated, leaving the spatial direction aligned with the vector (x = l,y = 1, z = 
1) distinct from any other directions. To sensibly describe an isotropic solution within the 
Bianchi I formulation, we should have no off-diagonal correlations; that is, instead of the 
formal isotropic condition (15.4)) . we should impose the physical isotropic condition: 



Vx = v 2 = v 3 = v, K 1 = K 2 = K 3 . 
G\i = G22 = G33 = G{2) but G12 
and similarly for other G^ n=2 \ 



K, 

- G13 



Si-- 
G23 



S2 

= o, 



S-t = s, 



(5.11) 



It is nontrivial to note that this prescription is also compatible with (14.381) — (14.491) . Particu- 
larly, the physical isotropic condition reduces (14.46)) to 



dv 1 
i— « ±37 v 



(K) 2 + (R) 2 G 



(2') 



(Giry^v 2 ' (67r 7 £ P1 ) 4 ^ 4 



1/2 



and 0S37D to 



Kl = QIC ( 1 - % \ . 



(5.12) 



(5.13) 



Consequently, we have [cf. (15.81) 



1 dp ^ , 16nG 
dpdcj) V 3 

with the critical value of given by [cf. ( 15.101) ] 



P±_ 

Pcrit . 



1/2 



(5.14) 



Pcrit 



1 ~ 4/3 
1 + 



Ppi < 3p P i. 



(5.15) 



The physical prescription (15. lip yields the effective equations slightly different from those 
by the formal prescription (15 ,4p . Comparing (15. 121) . (which gives isotropic solutions within 
the Bianchi I model) with (I5.5j) . (which describes the solution in the isotropic model), we 
note that, at least up to G^ n=2 \ the equation of motion for p subject to ( 15.111) receives no 
back-reaction due to the squeezing G$\, while that subject to (15.41) does. Furthermore, 
(15.10)) and (15 . 15)) do not agrees precisely. The moral we learned is: In the context of 
effective theory, it is possible to have a well-posed solution exhibiting a certain symmetry 
(e.g. isotropy) within the framework which does not assume such a symmetry; however, the 
well-posed solution in the less symmetric (e.g. Bianchi I) theory and the solution obtained 
directly from the more symmetric (e.g. isotropic) theory agree only approximately. 



VI. SUMMARY AND DISCUSSION 



When the effective theory developed in |13J is applied to the cosmological models more 
complicated than the k = FRW case, we face the problem that, in the scheme of Approach 



24 



I, the Hamiltonian operator is not polynomial of the fundamental operators and, as a con- 
sequence, the approximation scheme is out of good control. To work around the technical 
difficulties, we devise a new approach — Approach II — for the constrained quantum system 
that possesses an internal clock. In Approach II, the quantum evolution are approximated 
by a finite system of coupled but ordinary differential equations (12.111) and (I2.12p for the 
classical and quantum variables and, on top of them, the Hamiltonian constraint is weakly 
imposed by (12.131) . The evolution equations (12.111) and (12.121) adhered to (I2.13P can be 
understood as the natural extension of the classical Hamilton's equations and the classical 
Hamiltonian constraint. 

It is tantalizing that Approach II, albeit only heuristically motivated, could be used as a 
viable formulation for constrained quantum systems even in the fundamental level, as it is 
shown to be self-consistent. Approach II has the virtue that the technical and conceptual 
difficulties in the standard treatment for the constrained quantum systems are avoided and 
the quantum evolution is posed in a very intuitive picture. Furthermore, the philosophy of 
timeless formulation is better actualized, as the notion of measurement of time is retained 
and the time variable is on the equal footing as other observables, giving insight into the 
problem of time in the relativistic quantum mechanics. 

When applied to the WDW theory in the Bianchi I model, up to the 2nd order, Approach 
II gives the effective solutions (I3.32j) (I3.36j) subject to the constraints (13.471) and (13.481) . 
which agree with the results obtained in the fully developed quantum theory in [8J] for the 
case without off-diagonal squeezing. We also expect that the evolution in the fully developed 
quantum theory receives back-reaction arising from the off-diagonal squeezing. 

For simplified LQC, Approach II leads to the effective equations of motion (14 . 38 p (14 . 43 p 
subject to the constraints (I3.47j) and (13.481) . The big bounces take place up to three times, 
once in each direction of pi, whenever ( I4.46[) flips signs. The bouncing scenario predicted in 



the semiclassical approach in [10|, [12| is affirmed as the directional densities qi indicate the 



occurrence of bounces, but the critical values £>/ iCr it are modified by the quantum corrections 
as shown in (I4.56[) . The off-diagonal squeezing gives further back-reaction and makes the 
evolution of pi in different diagonal directions couple to one another. To see if these results 
agree with the standard treatment of the simplified LQC, it is necessary to explore more 
details of the numerical computation in [||. 

Additionally, the framework of effective equations offers a language to describe the 
isotropy-reduced solution in the context of the anisotropic Bianchi I model. The fact that 
the physical prescription for the isotropic condition is compatible with the effective equations 
of motions implies that, within the anisotropic framework, it is possible to have a well-posed 
solution exhibiting isotropy, which, however, is slightly different from the solution in the 
isotropic framework. The lesson may teach us how a certain symmetry emerges from a less 
symmetric model and whether a symmetric theory can be systematically derived from a 
more fundamental (namely, less symmetric) theory. 

However, we should keep the caveat in mind that the methodology of Approach II has not 
been stringently validated and it is still unclear how well and why Approach II agrees with 
the standard treatment for constrained quantum systems for generic cases. Therefore, a 
more rigorous formulation directly derived from the standard quantum theory is on demand 
and currently under development }18| . Meanwhile, to test robustness of ramifications of 
Approach II, it will be very instructive to repeat the analysis of this paper to other extended 
models, such as LQC in the k — ± FRW models (1, @], LQC in Kantowski-Sachs spacetime 
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[2~5| and the loop quantum geometry of the Schwarzschild black hole interior 26[. 14 

Finally, as commented in Sec. IIV Al the /i'-scheme quantization has a better scaling 
property but the quantum theory based on it is very difficult to construct. The difficulty 
is due to the fact that the corresponding affine variables do not exist. (See Appendix B of 
[l~2| for more details.) In the language of Approach II, the difficulty is translated into the 
problem that we cannot find a finite set of fundamental variables which form a closed algebra 
of Poisson brackets, as opposed to (I4.17p . (14.181) and (I4.19P for the fl scheme. Nevertheless, 
there might be a physical sense of approximation to truncate the infinite set of fundamental 
variables and then the treatment of Approach II can be carried over again. This strategy 
might gives a sound machinery to study the effective dynamics of the /i'-scheme LQC, even 
if the detailed construction of the quantum theory remains unknown. 
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APPENDIX A: THE WDW THEORY IN APPROACH I 

As a comparison to Approach II in Sec. IHIt the effective theory of Approach I for the 
WDW theory is presented here along the line of Sec. Ill Bl 

In Approach I, the Hilbert space of the WDW theory is 7Y WDW = L 2 (M. 3 ,d 3 p) and (f) is 
treated as a pure time variable. That is, |^ (</>)) G 7i WDW is a state at the instant cf) and the 
dynamics is govern by the Schrodinger equation 

-ih-l^\*(<j>)) = ±hV$\*(<j>)) (Al) 

as stated in (12.101) with the Hamiltonian operator given by 

± = + - 1 (K 2 K 3 + K 1 ^ + K l K 2 \ V2 (A2) 

A/47TG7 2 v ' 

as VG) is defined in (13.61) . 

Associated with the Hamiltonian operator in (1A2[) . the quantum Hamiltonian defined in 
(12T7D yields 

Hq = ±(n) (A3) 
with (p^) given by (I3.45j) in terms of classical and quantum variables. 



14 In the Schwarzschild black hole interior, the classical Hamiltonian constraint can be rescaled as H = 
— (2bcpbp c + (b 2 +7 2 )p^). Even though there is no matter content, pt could be locally treated as the 
internal clock by identifying the quadratic term 7 2 p 2 in H as pi. 
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Equation (12.51) gives the effective equations of motion: 



dK 1 

— — = {K 1 , H Q \ = K 1 := ^irG^hK, 1 are constant, (A4) 

dcp 

= {Pi,#q} ~T<7 Pi A5 



where * is short for 



i /tttt^ 1 , 3(if 2 + iY 3 ) 2 n (^)-3^)(^ + ir 3 ) 22 
x ^ 29(K 1 ) 1 / 2 16#(iY') 5 /2 16 0(^)5/2 (i^ + K 3 ) 

>2s 



OTzgg+gj r 33 

16 0(^)5/2(^2 + ^3) 

(if 2 + iY 3 )(iY l2 + K 2 K 3 ) + K\K 22 + 6iY 2 iY 3 + if 32 ) 23 
8 0(K 7 ) 5/2 (^ 2 + ^ 3 ) 



B{K') - 3iY 22 Gl3 _ OjK 1 ) - 3iY 32 Gl2 



1 + 



G 7 



2^2 



(A6) 



%6{K l Y/ 2 8 0(K r ) 5 / 2 

and so on for p7 2 , p^ 3 in the cyclic manner. Similarly, by ( 12.61) . we have 
dG 1J 

—— = {G IJ ,H Q } = G /J are constant, (A7) 



dG\\ 



{G 12 , H Q } = Tl ~iG 12 J^!±^! + &±&\ + 7 -i p -i (G(' 



(A8) 
(A9) 



^ = {G[,H Q } = Tl- 1 G[! [ ^^y^O(G^). (A10) 

Note that p^j are constants of motion and p^j ~ = (p^) if |G /J | <C SirGp^,. 

The resulting effective equations of motion are qualitatively different from those in (13.211) . 
(I3.22p and fj3.25j) — (13.28[) . Approach II agrees with Approach I only if all the 2nd order quan- 
tum variables are ignored. In particular, (I3.22|) receives back-reaction only if off-diagonal 
squeezing are present (Gj ^ for I ^ J) but (1A5j) receives back-reaction both from diagonal 
and off-diagonal squeezing. Furthermore, even with squeezing terms all vanishing, f]A5|) still 
receives back-reaction through p^j and gives a trajectory of (pi) slightly deviated from the 
classical one, which disagrees with the result of [s[ for the fully developed WDW theory. 

The comparison suggests that, for the Bianchi I model, the effective description of Ap- 
proach I is more sensible than that of Approach II. The clumsiness of Approach II is due 
to the non-polynomiality of (1A2|) . which spoils the order- by-order approximation scheme. If 
we formally reduce the Bianchi I model to the isotropic one by taking (15.41) . it can be easily 
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shown that p^i = p^ = p^z = (p^) and Approach II yields the same result as Approach 
I. That said, Approach I is good for the k = FRW model as the polynomiality of the 
Hamiltonian operator is recovered. 



APPENDIX B: SELF-CONSISTENCY OF APPROACH II 

In order to be a viable formulation for constrained quantum systems, Approach II has to 
pass two tests for self-consistency as mentioned Sec. Ill CI 

For the first test, to see if the weakly imposed constraint (12.131) still holds after the given 
initial state being evolved, we define the operators 

Q k := (±p 4> ) k - {h 2 Q) k ' 2 = (±pt) k - (-2H ) h / 2 (Bl) 

and the associated variables 

Q k := (Q k ) = (± P(t> ) k - (-2# ) fc/2 (B2) 

to see if Q k = remains. 
To check this, we compute 

= {n k ,H Q } = ±([h k ,H\) 

= ^({(±P<p) k -(- 2 Ho) k/2 ,Pl + 2H ]) 

= ±k ((±fc) {-2H Q ) k '^ V'm Wcyl =F k {{±ht~ l V'{$)) w , (B3) 

where V(<f>) G Ho represents the potential of cf). In the case with a free massless scalar field, 
V(<f) = and we have dfl k /dt' = exactly. If V(<p) ^ 0, the right-hand side of flB3p does 
not vanish except for the case of k = 2. Nevertheless, since we are only interested in the 
quantum evolution in the realm close to classical behavior, we can make the "factorization 
approximation" : 

(A 1 A 2 ...A n ) Weyl ^(A 1 )(A 2 )...(A n ) (B4) 

for hermitian operators as long as the quantum state is semiclassical enough. Consequently, 
(IB3I) can be approximated and gives 

^ « ±k{(±p,)((-2H ) k / 2 - 1 )(V>m - ((±^) fc - 1 )(^'(0))} « (B5) 

with the help of f|2.13p . Therefore, the first test is rigorously verified if V(<p) = and 
reasonably justified if V(4>) ^ 0. In particular, we can see why a free massless scalar field is 
special for being the internal clock. 

Next, for the second test, we consider the lapse function N = N(q l ,pi) associated with the 
time coordinate t' via dr = Ndt' with r being the proper time. Let H T be the (unsealed) 
Hamiltonian associated with r; the rescaled Hamiltonian associated with t! is given by 
H* = NH T . In quantum theory, the Hamiltonian operator corresponding to H l is IP = 
(NH T ) 

Weiy with N = N(q % ,Pi)w e \ y . With respect to t', the time derivative of an arbitrary 
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variable lo is given by 



« (iV)i(^,^]) + (^)i(^iV]) 

= (iV>{^}, (B6) 

where again the factorization approximation is made and we also use (H T ) = by (12.131) . 
In particular, we have 

^~(N){4>,HT Q } (B7) 

and subsequently the factor (N) cancels when we compute duo /dip. Therefore, the resulting 
effective dynamics, when related to the internal time 0, is independent of different choice 
of the lapse function . However, it is not guaranteed that the error due to the factorization 
approximation is always of higher order than the order of our interest. 



APPENDIX C: POISSON BRACKETS OF THE STANDARD VARIABLES 



If we use the canonical pairs c 1 := (c 1 ) and pi := (pi) as the classical variables, the 
associated quantum variables are defined as: 



G (n=2) . 



Gij 


■= ((Pi ' 


(vi))(vj - 


(PJ))), 


G IJ 


■= ((c 1 - 




<£ J »>, 


G l j 


:= <(c' - 


(&))(PJ- 


(Pj)))weyl 



(CI) 



This is the standard formalism, in which the quantum variables commute with the clas- 
sical variables. For reference, the Poisson brackets of the standard quantum variables of the 
2nd order are listed: 



{Gj, Gf} = (8vrG 7 ) 2 {SlG* - dfGi) , {G IJ , G KL ) = 4(8vrG 7 ) 2 ^G^, 
{G IJ ,G« } = 2(8vrG 7 ) 2 4 / G J)X , {G u ,Gf } = -2(87rG 7 ) 2 ^G J}i , 
{G LJ ,G KL } = {G IJ ,G KL } = 0. 
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